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Abstract. Earlier we proposed the stochastic point process model, which reproduces 
a variety of self-affine time series exhibiting power spectral density S{f) scaling as 
power of the frequency / and derived a stochastic differential equation with the same 
long range memory properties. Here we present a stochastic differential equation as a 
dynamical model of the observed memory in the financial time series. The continuous 
stochastic process reproduces the statistical properties of the trading activity and serves 
as a background model for the modeling waiting time, return and volatility. Empirically 
observed statistical properties: exponents of the power-law probability distributions 
and power spectral density of the long-range memory financial variables are reproduced 
with the same values of few model parameters. 
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1. Introduction 

Stochastic volatility models of the financial time series are fundamental to investment, 
option pricing and risk management ^1 The volatility serves as a quantitative 
price diffusion measure in widely accepted stochastic multiplicative process known as 
geometric Brownian motion (GBM). Extensive empirical data analysis of big price 
movements in the financial markets confirms the assumption that volatility itself is a 
stochastic variable or more generally the function of a stochastic variable |3] . By analogy 
with physics we can assume that speculative prices p{t) change in a "random medium" 
described by the random diffusion coefficient. Such analogy may be reversible from 
the point that the complex models of stochastic price movements can be applicable for 
the description of complex physical systems such as stochastic resonance, noise induced 
phase transitions and high energy physics applications. This analogy contributes to 
further development of statistical mechanics - nonextensive one and superstatistics have 
been introduced IHEj- 

Additive-multiplicative stochastic models of the financial mean-reverting processes 
provide rich spectrum of shapes for the probability distribution function (PDF) 
depending on the model parameters J)\. Such stochastic processes model the empirical 
PDF's of volatility, volume and price returns with success when the appropriate 
fitting parameters are selected. Nevertheless, there is the necessity to select the most 
appropriate stochastic models able to describe volatility as well as other variables under 
the dynamical aspects and the long range correlation aspects. There is empirical 
evidence that trading activity, trading volume, and volatility are stochastic variables 
with the long-range correlation [3 |H] and this key aspect is not accounted for in 
widespread models. Moreover, rather often there is evidence that the models proposed 
are characterized only by the short-range time memory [H]. 

Phenomenological descriptions of volatility, known as heteroscedasticity, have 
proven to be of extreme importance in the option price modeling Autoregressive 
conditional heteroscedasticity (ARCH) processes and more sophisticated structures 
GARCH are proposed as linear dependencies on previous values of squared returns 
and variances jTUl E]- These models based on empirically fitted parameters fail in 
reproducing power law behavior of the volatility autocorrelation function. We do 
believe that the stochastic models with the limited number of parameters and minimum 
stochastic variables are possible and would better refiect the market dynamics and its 
response to the external noise. 

Recently we investigated analytically and numerically the properties of stochastic 
multiplicative point processes [12] , derived formula for the power spectrum and related 
the model with the general form of multiplicative stochastic differential equation |13j . 
Preliminary the comparison of the model with the empirical data of spectrum and 
probability distribution of stock market trading activity jHj stimulated us to work 
on the definition of more detailed model. The extensive empirical analysis of the 
financial market data, supporting the idea that the long-range volatility correlations 
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arise from trading activity, provides valuable background for further development of the 
long-ranged memory stochastic models [ZllHl- We will present the stochastic model of 
trading activity with the long-range correlation and will investigate its' connection to 
the stochastic modeling of volatility and returns. 



2. Stochastic model of interevent time 



Earlier we proposed the stochastic point process model, which reproduced a variety of 
self-affine time series exhibiting the power spectral density S{f) scaling as power of the 
frequency / P^ll4j . The time interval between point events in this model fluctuates as 
a stochastic variable described by the multiplicative iteration equation 

Tk+i = rk + irl^'^ + ar^Ek- (1) 

Here interevent time = tk+i—tk between subsequent events k and k + 1 fluctuates due 
to the random perturbation by a sequence of uncorrelated normally distributed random 
variable {sk} with the zero expectation and unit variance, a denotes the standard 
deviation of the white noise and 7 -C 1 is a coefficient of the nonlinear damping. It has 
been shown analytically and numerically ^3 ^^^^ point process with stochastic 
interevent time may generate signals with the power-law distributions of the signal 
intensity and 1//^ noise. The corresponding Ito stochastic differential equation for the 
variable r(t) as a function of the actual time can be written as 

dr = -fT^^'-^dt + ar^-'/^dW, (2) 

where is a standard random Wiener process. Eq. (0) describes the continuous 
stochastic variable r(t) which can be assumed as slowly diffusing mean interevent time 
of Poisson process with the stochastic rate l/r(t). We put the modulated Poisson 
process into the background of the long-range memory point process model. 

The diffusion of r must be restricted at least from the side of high values. Therefore 
we introduce a new term — ycr^ (^^^ -j-'^i^-'^ into the Eq. which produces the 

exponential diffusion reversion in equation 



dr 



m 2 
7 — 77'^ 



T^^'-Mt + (XT^-^/^dW, 



(3) 



where m and Tq are the power and value of the diffusion reversion, respectively. The 
associated Fokker-Plank equation with the zero flow will give the simple stationary PDF 

T 



Pr 



~ r°+^ exp 



(4) 



with a = 2(70-— /i), where 7o- =7/(1^. We define the conditional probability of interevent 
time Tp in the modulated Poisson point process with stochastic rate 1/r as 

(5) 



1 ^ ^ 




rj = — exp 
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Then the long time distribution f{Tp) of interevent time Xp has the integral form 



(^(rp) = C 
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with C defined from tlie normalization, <^(rp)drp = 1. In tlie case of pure exponential 
diffusion reversion, m = 1, PDF © has a simple form 

^ir,) = I (^) K(i+.) (2 ^] , (7) 
r(2 + a)ro \ToJ V V ^0 / 

where [z) denotes the modified Bessel function of the second kind. For m > 1 
more complicated structures of distribution f{Tp) expressed in terms of hypergeometric 
functions arise. 



3. Stochastic model of flow of points or events 

The introduced stochastic multiplicative model of interevent time, the interval between 
trades in the financial market, defines the model of event flow n. First of all we apply 
Ito transformation of variables introducing flow of events nit) = The stochastic 
differential equation for n follows from Eq. ((21) , 

n^^'^dt + an'^dW, (8) 

where r] = ^ — fi and uq = I/tq. Eq. (jS)) describes stochastic process n with PDF 

P(n)~i^exp{-(^)"}, A = 2(r/-l + 7.), (9) 

and power spectrum S{f) jT^llISllIll 

Noteworthy, that in the proposed model only two parameters, 70- and rj (or fi), define 
exponents A and (3 of two power-law statistics, i.e. of PDF and power spectrum. Time 
scaling parameter cx^ in Eq. (jSI) can be omitted adjusting the time scale. 

Stochastic variable n denotes the number of events per unit time interval. One has 
to integrate the stochastic signal Eq. (jH)) in the time interval to get number of events 
in the selected time window. In this paper we will denote the integrated number of 
points or events as 

N{t,T^) = J'^^\{t')dt' (11) 

and will call it trading activity in the case of the financial market. Flow of points 
or events arises in different fields, such as physics, economics, cosmology, ecology, 
neurology, the Internet, seismology, i.e., electrons, photons, cars, pulses, events, and so 
on, or subsequent actions, like seismic events, neural action potentials, transactions in 
the financial markets, human heart beats, biological ion-channel openings, burst errors 
in many communication systems, the Internet network packets, etc. We will discuss 
possible application of the proposed stochastic model to model the trading activity in 
the financial markets. 
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4. Stochastic model of trading activity 

It is widely accepted that in high-frequency financial data not only the returns but also 
the waiting times between the consecutive trades are random variables ^S]. Waiting 
times between trades do not follow the exponential distribution and the related point 
process is not the Poisson one. The extensive empirical analysis provides evidence that 
the related stochastic variable trading activity defined as fiow of trades is stochastic 
variable with the long range memory PS]- We will investigate how the proposed 
modulated Poisson stochastic point process can be adjusted to model trading activity 
with the empirically defined statistical properties. Detrended fiuctuation analysis is 
one of the methods to define the second order statistics, the autocorrelation of trading 
activity. The histogram of the detrended fiuctuation analysis exponents z/ obtained by 
fits for each of the 1000 US stocks shows a relatively narrow spread of i' around the 
mean value u = 0.85 ± 0.01 JH]- We use relation between the exponents of detrended 
fiuctuation analysis and the exponents of power spectrum /3 = 2z/ — 1 ^\ and in this way 
define the empirical value of the exponent for the power spectral density (3 = 0.7. Our 
analysis of the Lithuanian stock exchange data confirmed that the power spectrum of 
trading activity is the same for various liquid stocks even for the emerging markets |18j . 
The histogram of exponents obtained by fits to the cumulative distributions of trading 
activites of 1000 US stocks [_l6j gives the value of exponent A = 4.4 ± 0.05 describing 
the power-law behavior of the trading activity. Empirical values of /5 = 0.7 and A = 4.4 
confirm that the time series of the trading activity in real markets are fractal with the 
power law statistics. Time series generated by stochastic process (jH)) are fractal in the 
same sense. 

Nevertheless, we face serious complications trying to adjust model parameters to 
the empirical data of the financial markets. For the pure multiplicative model, when 
/i = 1 or ?7 = 3/2, we have to take 7o- = 0.85 to get P = 0.7 and 7o- = 1.7 to get 
A = 4.4, i.e. it is impossible to reproduce the empirical PDF and power spectrum with 
the same relaxation parameter and exponent of multiplicativity /z. We have proposed 
possible solution of this problem in our previous publications ^1 1^ deriving PDF for 
the trading activity 



When 7 ^ this yields exactly the required value of A = 5 -|- 2a = 4.4 and 

p = 0.7 for 7^ = 0.85. 

Nevertheless, we cannot accept this as the sufiiciently accurate model of the 
trading activity as the empirical power law distribution is achieved only for very high 
values of the trading activity. Probably this reveals the mechanism how the power 
law distribution converges to normal distribution through the growing values of the 
exponent, but empirically observed power law distribution in wide area of N values 
cannot be reproduced. Let us notice here that the desirable power law distribution of 
the trading activity with the exponent A = 4.4 may be generated by the model (jH)) 
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Figure 1. Power spectral density S{f) calculated with parameters 7 — 0.0004; 
a = 0.025; e = 0.07; tq = 1; m = 6. Straight lines approximate power spectrum 
S ^ with Pi = 0.33 and (32 = 0.72: a) S{f) calculated by the Fast Fourier 

Transform of n series generated by Eq. H13|l . b) S{f) averaged over 20 series of 100000 
iterations of the flow I{t) — J2k ^i^ " ^k) with the interevent time — t^+i — tk 
generated by Eq. (|14|l . 



with f] = 5/2 and 7o- = 0.7. Moreover, only the smallest values of r or high values 
of n contribute to the power spectral density of trading activity |TB]. This suggests 
us to combine the point process with two values of fi: (i) /i ~ for the main area 
of diffusing r and n and (ii) /i = 1 for the lowest values of r or highest values of 
n. Therefore, we introduce a new stochastic differential equation for n combining two 
powers of multiplicative noise. 



where a new parameter e defines crossover between two areas of n diffusion. The 
corresponding iterative equation of form for in such a case is 



Eqs. (IT!^ and (fT^ define related stochastic variables n = ^ and r, respectively, and 
they should reproduce the long range statistical properties of the trading activity and 
of waiting time in the financial markets. We verify this by the numerical calculations. 
In figure HI we present the power spectral density calculated for the equivalent processes 
(fT^ and (fT^ (see [14j for details of calculations). This approach reveals the structure 
of the power spectral density in wide range of frequencies and shows that the model 
exhibits not one but rather two separate power laws with the exponents f3i = 0.33 and 
(32 = 0.72. From many numerical calculations performed with the multiplicative point 
processes we can conclude that combination of two power laws of spectral density arise 
only when multiplicative noise is a crossover of two power laws, see (|T3|l and (|Tl|l . We 
will show in the next section that this may serve as an explanation of two exponents of 
the power spectrum in the empirical data of volatility for S&P 500 companies ^H] • 

Empirical data of the trading activity statistics must be modeled by the integrated 
flow of event N defined in the time interval ^ tq. In figure El we demonstrate the 
cumulative probability distribution functions P>{n) calculated from the histogram of 
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Figure 2. Cumulative PDF P>{n) calculated from the histogram of N/t^ generated 
by Eq. H14() with the increasing time interval Td from the above curve: Td = 1, 10, 50 
and 250. Other parameters are as in figure ^ 



N/Td generated by Eq. (fT^ with increasing time interval ra. This illustrates how 
distribution of the integrated signal converges to the normal distribution (the central 
limit theorem) through growing exponent of the power-law distribution and provides an 
evidence that the empirically observed exponent A = 4.4 of the power-law distribution 
of in IH] can be explained by the proposed model with the same parameters suitable 
for description of the power spectrum of the trading activity. 

The power spectrum of the trading activity N can be calculated by the Fast Fourier 
Transform of the generated numerical series. As illustrated in figure IHl the exponents 
P = 0.7 of the power spectrum are independent of and reproduce the empirical results 
of the detrended fluctuation analysis [3 |H] . 

The same numerical results can be reproduced by continuous stochastic differential 
equation (fT^ or iteration equation (fT^ . One can consider the discrete iterative equation 
for the interevent time Tk fll4|l as a method to solve numerically continuous equation 
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(15) 



The continuous equation (fT^ follows from the Eq. (fT3j) after change of variables 
n = 1/t. 

We can conclude that the long range memory properties of the trading activity in 
the financial markets as well as the PDF can be modeled by the continuous stochastic 
differential equation (fT^ . In this model the exponents of the power spectral density, j3, 
and of PDF, A, are defined by one parameter 70- = 7/0"^. We consider the continuous 
equation of the mean interevent time r as a model of slowly varying stochastic rate 1/t 
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Figure 3. Power spectral density of the trading activity N calculated by the Fast 
Fourier Transform of N series generated with Eq. (|13|l for the same parameters as in 
figures Hand 121 a) Td = 10; b) Td = 50; c) Td = 250. Straight lines approximate power 
spectrum S ^ 1//^, with /3 = 0.7. 
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Figure 4. Probability distribution function P{Tp) calculated from the histogram of 
Tp generated by Eq. (|16|l with rate calculated from Eq. H15(l . Used parameters are 
7 = 0.0004; a = 0.025; e = 0.07; tq = 1 and m = 6. Dashed line approximates power 
lawP(rp)^rp-0-i5. 



in the modulated Poisson process 
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In figure m we demonstrate the probabihty distribution functions -P(Tp) calculated 
from the histogram of Tp generated by Eq. with the diffusing mean interevent time 
calculated from Eq. (fT3j). 
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Numerical results show good qualitative agreement with the empirical data of 
interevent time probability distribution measured from few years series of U.S. stock 
data This enables us to conclude that the proposed stochastic model captures the 
main statistical properties including PDF and the long range correlation of the trading 
activity in the financial markets. Furthermore, in the next section we will show that 
this may serve as a background statistical model responsible for the statistics of return 
volatility in widely accepted GBM of the financial asset prices. 

5. Modeling returns and volatility 

We follow an approach developed in jTj |H1 E] to analyze the empirical data of price 
fluctuations driven by the market activity. The basic quantities studied for the individual 
stocks are price p{t) and return 

xit,Td) = \np{t + Td)-\np{t) (17) 

Return x{t, ra) over a time interval can be expressed through the subsequent 
changes 6xi due to the trades ^ = 1, 2....N{t, tj) in the time interval [t, t + t^], 

xit,Td)= Yl (18) 
1=1 

We denote the variance of 6xi calculated over the time interval as W'^(t,T^). 
If 6xi are mutually independent one can apply the central limit theorem to sum (fTH|l . 
This implies that for the fixed variance W^{t, Td) return x(t, t^) is a normally distributed 
random variable with the variance W^{t, Td)N{t,Td) 

x{t, Td) = W{t, Td)^jN{t,Td)et. (19) 

Empirical test of conditional probability P{x{t,Td)\W(t,Td)) [T confirms its' 
Gaussian form, and the unconditional distribution P(x(t, Td)) is a power-law with the 
cumulative exponent 3. This implies that the power-law tails of returns are largely due to 
those oiW{t,Td)- Here we refer to the theory of price diffusion as a mechanistic random 
process 122] • For this idealized model the short term price diffusion depends on 
the limit order removal and this way is related to the market order flow. Furthermore, 
the empirical analysis confirms that the volatility calculated for the fixed number of 
transactions has the long memory properties as well and it is correlated with real time 
volatility jSHj. We accumulate all these results into strong assumption that standard 
deviation W{t,Td) may be proportional to the square root of the trading activity, i.e., 
W(t, Td) = k^N{t, Td). This enables us to propose a very simple model of return 

x(t, Td) = kN{t, Td)et (20) 

and related model of volatility v = |a;(t, rd)| based on the proposed model of trading 
activity (fT^ . We generate series of trade flow n{t) numerically solving Eq. ()13p 
with variable steps of time Ati = hi = no/rii and calculate the trading activity in 
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Figure 5. (a) Cumulative probability distribution function of the volatility, P>(w), 
averaged over 10 intervals calculated from the series of n{t) generated by Eqs. ifT^ and 
()20|l . all parameters are the same as in previous calculations. Dashed line approximates 
the power law P{v) ~ ^/v'^^^. (b) Power spectral density S{f) of w calculated from FFT 
of the same series n{t). Straight lines approximate power spectral density S ^ 1/ f^^-^ 
with (3i = 0.6 and (32 = 0.24. 

subsequent time intervals as N{t, Td) = //^^"^ n{t')dt'. This enables us to generate 
series of return Td), of volatility f(t, Td) = |x(t, rd)| and of the averaged volatility 

In figure El we demonstrate cumulative distribution of v and the power spectral 
density of v(t, t^) calculated from FFT. We see that proposed model enables us to catch 
up the main features of the volatility: the power law distribution with exponent 2.8 and 
power spectral density with two exponents /?i = 0.6 and /?2 = 0.24. This is in a good 
agreement with the empirical data [T^ I23j. 

6. Conclusions 

Earlier proposed stochastic point process model |14[ ITH] as a possible model of trading 
activity in the financial markets has to be elaborated. First of all, we define that 
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the long-range memory fluctuations of trading activity in financial markets may be 
considered as background stochastic process responsible for the fractal properties of 
other financial variables. Waiting time in the sequence of trades more likely is double 
stochastic process, i.e., Poisson process with the stochastic rate defined as a stand-alone 
stochastic variable. We consider the stochastic rate as continuous one and model it 
by the stochastic differential equation, exhibiting long-range memory properties. We 
reconsider previous stochastic point process as continuous process and propose the 
related nonlinear stochastic differential equation with the same statistical properties 
[T^ . One more elaboration of the model is needed to build up the stochastic process 
with the statistical properties similar to the empirically defined properties of trading 
activity in the financial markets. We combine the market response function to the noise 
as consisting of two different powers: one responsible for the probability distribution 
function and another responsible for the power spectral density. The proposed new 
form of the continuous stochastic differential equation enables us to reproduce the 
main statistical properties of the trading activity and waiting time, observed in the 
financial markets. More precise model definition enables us to reproduce power spectral 
density with two different scaling exponents. This provides an evidence that the market 
behavior is dependant on the level of activity and two stages: calm and excited must 
be considered. We proposed a very simple model to reproduce the statistical properties 
of return and volatility. More sophisticated approach has to be elaborated to account 
for the leverage effect and other specific features of the market. 
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